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SECTION I 
SUMMARY 


Objectives 

The objectives of the first contract "Study of the Liquid-Solid Tran- 
sition for Materials Processing in Space" have been discussed in 
Reference 1, and these objectives continue to apply to the second year's 
effort (Modification 1) which is described in this report. The objectives 
are: 


1. To analyze the behavior of dense liquids near the 
solidification point while the liquid in question is 
under the influence of magnetic fields or near-zero 
gravity conditions, and 

2. to do this within the framework of existing liquid 
state models and classical field theory. 

In the present work, these objectives have been extended to include 
numerical calculation of the magnitudes of external field effects on 
liquids and to describe how external fields can influence the sub- 
structure of the solid during the solidification process. The materials 
considered in this analysis are (as much as possible) semiconducto s 
and metals, since these are the materials of most interest in the Space 
Processing Program. 

Approach 

The approach employed for this second year of the Liquid-Solid Transition 
Study was to build on the knowledge of liquid state models and field- 
liquid interactions gained during the first year. This was done by 
defining three tasks and several subtasks: 

Task 1 - Liquid Property Data Search 
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Task 2 

2A 

2B 

2C 


Determination of Field Effects on Solute 
Distri bution 

Analysis of Field Distributions in Melts 

Calculation of Field Forces on Solute Atoms 
and Resultant Velocities 

Relation of Microsegregation to Solute Atom 
Velocity 


Task 3 
3A 
3B 

3C 


Analysis of Field Effects on Liquids 

Calculation of Magnetoviscous Effects 

Calculation of Field Effects on Diffusion and 
Related Phenomena 

Analysis of Body Forces Due to Oscillating 
Magnetic Fields 


The experimental values of basic liquid properties obtained in Task 1 
were used in Task 3 «. input to the theoretical equations derived in 
Reference 1. Calculations performed in T*sk 3 then provide quantitative 
estimates of the effects which magnetic and gravitational fields have 
on melts of metals or semiconductors. Task 2 was essentially a new 
effort to develop a model of the solidification process which would 
describe to some extent the effects external fields will produce on 
the substructure of the solid. 


The qualitative analysis performed during the initial year of the 
study suggested approaches for further investigation and calculation 
in several areas. The magnetodynamic effects of interest included the 
magnitude of the induced magnetic viscosity and its importance relative 
to the molecular viscosity as measured by the Hartmann number. Analysis 
of radio-frequency (R-F) field effects can easily be generalized to 
any oscillating magnetic field since all time-varying magnetic fields 
induce eddy currents in conducting media. These eddy currents then 
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interact with the magnetic field to produce body forces on melts of 
metals or semiconductors placed in an oscillating magnetic field. In 
order to obtain a quantitative feeling for these eddy current body 
forces, they are compand to both the viscous force induced by the 
field and to the gravitational body force, pg, as described in 
Section 2.2. External field effects on diffusive phenomena were 
approached through the Free Volume Model descriptions of diffusion 
coefficient and solidification rate. This approach describes only the 
direct effects fields have on these and related quantities, and neglects 
any effects of convection which are considered to some extent in Section 
2.2. The Task 2 description of how external fields affect solute 
segregation during solidification was approached through a four step 
procedure which (1) relates the field distribution inside a melt to 
the external field inducing it, (2) calculates the force on a solute 
atom due to the internal field, (3) relates solute atom velocity in 
the melt to the driving force found in step (2), and (4) uses Sekr ka ' s 
stability theory to show how solute atom velocity in the melt influences 
solute microsegregation in the solid. Thus emphasis was .aced on 
development of procedures for doing practical calculations of solidifi- 
cation parameters affected by external fields (in Task 3) and on 
development of a mathematical model which formally describes field 
effects on microsegregation (in Task 2). 

Results 

The search for experimental values of liquid property data has identified 
39 papers (Table I) and 12 books (Table II) containing data useful in 
the study of external field effects on melts of semiconductors or 
metals. Additional references dealing with liquid state theory and 
problems are given in Section 2.5, and each of these references them- 
selves contain references to works dealing with other problems or data 
pertinent to the study of real liquids and liquid state models. The 
results of the data search which apply tc the investigation of field 
effects on solidification are condensed in Tables III and V which 
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contain the input data for calculation of magnetic and gravitational 
field effects on diffusion coefficient, solidification rate and for 
calculation of field forces on individual molecules in the melt. 

There are two reasons why data for more than four materials does not 
appear in Table III. First, there has been very little experimental 
work done on liquid metals and semiconductors as compared to solids, 
gases and organic liquids, and so there is a dearth of experimental 
data for the physical properties of these materials. Second, for 
a given material, one must find some eleven material properties just 
to do the calculations required for the diffusion coefficients and 
solidification rates. The materials listed in Table III were the 
only liquids for which values of all eleven properties could be found 
simultaneously. 

Calculations of magnetic viscosity and Hartmann number show that 
magnetovi scous effects predominate in the metals and semiconductors 
studied for fields on the order of about 10 gauss or above. This is 
in agreement with the finding that gravity-driven convection velocity 
decays exponentially with a time constant inversely proportional to 
Me square of the magnetic induction field applied to the melt to 
induce the magnetic viscosity. Eddy current body forces induced in 
melts by oscillating magnetic fields are quadratic functions of 
field strength. Because the susceptibilities of most paramagnetic 
and diamagnetic liquids are on the order of 1(T 6 cgs units, the 
valuer, of the eddy current body force (f g ) for any of these liquids 
at a fixed field strength are equal within 0.001%. Comparison of f 
to other forces on the melt shows that eddy current forces usually 
predominate, even over the magnetovi scous force, which arises from 
the same field as f . This predominance depends, however, on the 
value of u, the initial convection velocity, in the magnetic Vis- 

142 

cosity case, and on the ratio — in the comparison of f g to 
gravitational body forces, where H is magnetic field strength and 
o is liquid density. 
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Turnbull's Free Volume Model of the liquid state and the liquid-solid 
transition was used to calculate direct effects of external fields 
on diffusion coefficients and solidification rates of selected 
representative materials. Direct effects are those derivable directly 
from the fields and not attributable to convective or magnetohydro- 
dynamic effects. The results of this analysis indicate that micro- 
gravity causes infinitesimal changes ( 10~ 5 % - 1 0“ 6 %) in diffusion 
coefficient and solidification rates whi^e the changes caused by a magnetic 

field of 10 5 oersteds are small (0.3% - 1.9%) but finite. Changes 
in diffusion coefficient or solidification rate in a magnetic field 
are dependent on the square of the field strength, and thus vary 
widely with field strength. The diffusion coefficient was found 
to increase in microgravity over its value on Earth, while whether 
it increases or decreases in a magnetic field depends on several 
parameters, primarily the field strength. The free energy term 
(see equation 14) decreases for all materials considered in micro- 
gravity, and its change in a magnetic field is dependent on the sign 
of the susceptibility difference, * x $ > and on the ratio of the 
free energy in zero field to the heat of fusion. Changes in solidifi- 
cation rate are, from equation (14) , dependent on the relative magnitudes 
of the free energy term change and diffusion coefficient change. In 
the examples chosen, the solidification rate decreased in microgravity 
and increased in a magnetic field of 10 s oersteds. 

It should be noted here that Wang computer programs have been written 
for calculation of the following parameters: 

Eddy current body force 

Ratio of eddy current force tc gravitational body force 
D'/D (magnetic field case) 

Free energy term (magnetic field case) 

Isothermal compressibi lity 
Coefficient of thermal expansion 
Entropy 

Specific heat at constant pressure 
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Free Volume Model a 2 parameter 

Conversion of density to specific volume. 

In Task 2 of the present contract, a model has been developed which 
formally relates external fields to solute distribution in solids. 

The existence of such a formal relationship constitutes a significant 
result, but the model in its present form is severely limited in its 
range of applicability and accuracy. The major limitations are as 
follows: 

(a) Exact, analytical solutions for magnetic field distri- 
butions satisfying the boundary conditions on melts 
for Czochralski crystal growth geometry cannot be 
obtained, 

(b) Random acceleration (gravity) fields in orbiting 
spacecraft are not well known, 

(c) Intermolecular forces in liquids are not well understood, 
and the problem of describing the effects of more than 
one field acting simultaneously on a molecule in a 
liquid requires more study, 

(d) Sekerka's theory of interface stability is capable of 
estimating segregation cell size to only an order of 
magnitude because the theory is based on the small 
perturbation approximation, and 

(e) Analysis of cases in which solute molecule velocity 

is in any direction other than that of the solidification 
velocity is extremely difficult and will require much 
more research. 

Results obtained under the above limitations include the determination 
of the occurrence of microsegregation as a function of the parameters 
V - j- where V is macroscopic solidification velocity, u is solute 
molecule velocity and p * 1-k with k being the segregation coefficient 
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for the material considered, and estimates of the size of cells of 
different solute concentration values in the solid. 

During the years this contract has been active, work performed under 
the contract has resulted in seven scientific publications (two 
already published i.i international jo: rnals) and four presentations 
to scientific meetings (one invited). These are listed in Appendix A. 

Conclusions and Recommendations 

The major conclusions of this study are: 

1. Values of the physical properties of molten semiconductors are 
difficult to obtain from the literature because very little 
experimental research has been done in this area, 

2. Of the research on liquid semiconductor properties which has 
been reported, much of it has been done by Russian researchers, 

3. Moderate and high magnetic fields produce strong forces in 
melts and damp out convective flows exponentially, 

4. Oscillating magnetic fields induce tody forces even stronger than 
the static field magnetoviscous force, 

5. The direct effects of microgravity on melts of the materials 
considered in this study were infinitesimal, while the effects 
of magnetic fields were dependent on field strength and were 
small but finite, 

6. The lack of a large direct effect of microgravity on liquids 
supports the thesis that the primary benefits of materials pro- 
cessing in space arise from secondary effects, such as the 
suppression of convection or the reduction of contamination by 
containerless processing. 
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7. A formal relationship between external fields and solute distribution 
in a solid during solidification has been established, 

8. The occurrence of microsegregation depends on the parameters 

(V - — ) and w, where u is related to the interface perturbation 
or cell size, 

9. There are several practical limitations to the theory which describes 
microsegregation dependence on external fields, including 

- the difficulty of obtaining gravitational (acceleration) field 
distributions in melts aboard an orbital laboratory or magnetic 
f ield distributions in melts during Czochralski growth on Earth, 

- the difficulty of representing the process of several fields 
acting on a molecule at once, and 

- the failure of the small perturbation approximation (which 
linearizes the temperature and concentration field equations) 
to yield results which agree well with experiment. 

Thus it is seen that, for the materials for which experimental property 
data are available, the similarity of microgravity and magnetic field 
effects on melts holds only for very low magnetic fields, and the 
effects are quantitatively quite different at higher magnetic field 
strengths. This was to be expected, since the magnetic interaction 
is much stronger than gravity, and has now been verified. 

The relationship of microsegregation in solids to external fields 
applied during solidification is found to be very complicated and 
not well understood at present. It is reconmended that much more 
research be done in this area, specifically in 

- representing field distributions in melts, 

- defining resulting forces on individual molecules and 
deriving the corresponding velocities, and 

- overcoming the limitations of the small perturbation 
approximation in calculating segregation cell sizes when 
microsegregation occurs. 
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It is further recommended that research begin on a flight experiment 
whose objective would be to compare the theory of Section 2.4 to 
experimental measurements of interface stability in a convection- 
free environment where the small perturbation approximation may be 
valid. 
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SECTION 2 

ANALYTICAL RESULTS 

2.1 LIQUID STATE PROPERTY DATA SEARCH 

Tne work detail eo in Reference 1 provides a qualitative theoretical 
description of many effects which magnetic and gravi .ational fields 
produce on melts during solidification processes, but there are two 
reasons why it is desirable to obtain a quantitative understanding 
of these effects. First, in most cases, the equations describing the 
field effects of interest do not give ar. obvious indication of the 
magnitude or direction of the change in liquid properties due to the 
field change without actually performing the calculations defined by 
the equations. Second, it is always helpful in understanding physical 
phenomena described by equations to become familiar with the numerical 
value of the various parameters in the equations as they correspond 
to different materials or conditions. For these reasons it was decided 
to perform numerical calculations of the field effects derived in 
Reference 1, and to do this, it was necessary to obtain numerical 
values of all parameters appearing in the equations describing the 
effects. 

These experimental values were to be obtained for materials similar 
to those most prominent in the Space Processing Program, i.e., those 
which show promise of improvement in tneir economically beneficial 
properties when processed in space. Semiconductors and metals, and 
in particular, Group III - Group V compounds were the materials 
receiving the most attention in tne Space Processing discipline at 
the beginning of the contract, thus the search for experimental data 
centered on these types cf materials. The search encompassed a com- 
puter search of NASA and DOL data banks, searches of the Physics 
Abstracts, the Science Citation Index, the Redstone Scientific Infor- 
mation Center book collection and personal contacts with experimental 
researchers at universities and other companies. It was found that 
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2.1 (Continued) 

much of the property data for molten semiconductors has been measured 
by Russian researchers, and that most of the Russian work deals with 
electrical properties of melts other than III-V compounds. Data for 
the III-V compounds corresponding to parameters and properties in the 
field effect equations^ proved extremely difficult, and in many 
cases impossible, to find. Tables I and II list references for the 
journal articles and books containing liquid property data which 
were identified in the search. The tables also indicate the particular 
materials and properties considered by each reference. It must be 
noted that each reference shown in these tables also contains references 
to other works on liquid metals or semiconductors, therefore these 
lists of references are intended as a starting point for researchers 
interested in a variety of materials and properties. Sufficient data 
was found for indium antimonide (InSb), germanium (Ge), mercury (Hg) 
and sodium (Na) to allow calculation of external field effects on 
temperature, volume, Gibbs free energy, diffusion coefficient and 
solidification rate as will be described in Section 2.3 (see Table III). 

2.2 MAGNETODYNAMIC EFFECTS uN MELTS 

Several macroscopic dynamic effects which are produced by magnetic 
fields are worthy of understanding in their own right and also may be 
helpful in understanding microgravity effects on melts. 

Calculations of Magnetic Viscosity and Hartmann Number 

Previous work has shown ^ that the magnetic viscosity induced in a melt 
with electrical conductivi ,j a is related to the magnetic induction $ 
by 

n m = oB ? -L 2 (1) 
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2.2 (Continued) TABLE I 

LIQUID STATE PHYSICAL PROPERTIES REFE^NCE LIST 



REFERENCE 

DATE 


PROPERTIES 

1. 

Glazov & Vertman 

1958 

InSb, GaSb. AlSb 

X 

2. 

Glazov A Petrov (in Russian) 

1958 

InSb, GaSb, AlSt 

rt 

3. 

Amirkhanov A Magomedov 

1964 

InSb 

r. 

4. 

Glazov A Chizhevskaya 

1963 

Ge, Si, ZnS 

X 

5. 

Blum A Ryabtsova 

1958 

InSb, GaSb 

0 

6. 

Busch A Vogt (in German) 

1954 

InSb 

a 

7. 

Glazov A ChizhevsKaya 

1961 

Ge, AlSb, GaSb, .iSb 

a 

8. 

Glazov A Chizhevskaya 

1962 

GaAs, InAs 

a i n 

9. 

Regel et.al. 

1971 

Sb 2 Se 3 

o ft R 

10 . 

Mogilevskii et.al. 

1971 

Sb 2 Te 3" Sb 2 Se 3 

a ft k 

11. 

Mogilevskii et.al. 

1972 

Sb, Te, Bi A Pb alloys 

thennodlf fusion 

12. 

Mogi evskii et.al. 

1972 

Tfe X Se l-X 

0 I K 

13. 

Fedorov A Machuev 

1971 

T1S, T1 2 S 

K 

14. 

Fedorov A Machuev 

1971 

GaSe, InSe 

K 

15. 

Mavlonov 

1972 

InSe, Ga 2 Sej 

o A k 

16. 

Kazandzhan et.al. 

1972 

Te-TI alloys 

o A TEP 

17. 

Kazandzhan et.al. 

1971 

I: A Se Impurities In 
T1 2 Te 

o A TEP 

18. 

Hodgkinson 

1970 

In-S A Tl. . 

o A theory 

19. 

Krestovnikov et.al . 

1967 

nost semiconductors 

av at melting 

20. 

Mai ' agov A Magomedov 

1969 

CuSbSe 2 

o, k A TEP 

21. 

Fedorov A Stll'bans 

1967 

Bi , Te, Sb, Se A Cu 9 
TeS L 

1C 

22. 

Fedorov et.al. 

1968 

ll 2 S« 3 

0, k A TEP 
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2.2 (Continued) 


TABLE I (Continued) 


REFERENCE 

DATE 

23. Fedorov A Machuev 

1967 

24. Fedorov & Machuev 

1968 

25. Glazov et.al. 

1967 

26. Perron 

1970 

27. Cutler & Mai Ion 

’964 

2E. Enderby & Walsh 

1966 

29. Enderby & Simnons 

19o9 

30. Cutler & Mai Ion 

1966 

31 . Edmond 

IS ' 

32. Ioffe A Regf ' 

1960 


33. 

Regel. et.al 

1970 

34. 

Regel , Smirnov. 

1972 


Shadrichev 


V, 

ison A Readall 

1962 


MATERIALS 

PROPERTIES 

Sb - Se alloys 

K 

SbS 

K 

AlSb. GaSb. InSb. InAs, 
GaAs 

Tm, H^, 

Te 1-x Se x 

O, it 

T1 x Te 1-x 

o A T £P 

CdSb, ZnSb, B1,Te„ 
Sb 2 Te 3 2 3 

R, TEP A o 

AupTe, CuTe, AgTe, Bl-Te, 
Tl-Te 

R, o 

Tl-Te solutions 

o, Seebeck coef 

As,Se, — As»Te, A 

o, TEP, optical 

As 2 Se 3 * T1 2 Te 3 

absorption 

NaCI . k. . . T1 S 2 , ^2^3* 
Sb 2 S 3> V 2 0j, Cu 2 S-FeS, 

SI. Ge.GaSb, InS\ HgTe. 
HgSe. CdTe, Se, i ,TeSe, 
BigTe CuTe, Cc ,5e, ZnSb 

*>* P» R* n 

TEP. mobility, i 

Tl-Se, Tl-Te. T1 2 Te, Tl 2 Se, 
T1 2 S. Sb 2 Se ^ , CuSbSe 2 

o , R. TEP 

Sb 2 Se 3 , Ge, SI, InSb. GaSb, 
ZnSb, B1 2 Te- , GeTe, SnTe, Te, 

o, k, L/L fl 

Cu 2 SnTe 3 , PbTe, Cu 2 GeTe 3 , GaTe 
T1 x Te y , CuT1Te 2 

9 

U’ 2 S 

1C 
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2.2 (Continued) 

TABLE 



REFERENCE 

PATE 

36. 

Johnson & Readal 1 

1963 

37. 

jtler & Hal Ion 

1962 

33. 

Cutler 6 Hal Ion 

1965 

39. 

Llchter 1 Soomelet 

1969 
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HATERIALS 

PROPERTIES 

Cu 2 S .25 Te .75* Cu 2 S 

K 

Te, Te-Se alloys 

k, o, TEP 

Tl-Te alloys 

TEP. P e 

InSb, GaSb, AlSb, InAs, 
GaAs 

' ”298* 
' S 298* 
6 T' H °298 , 


T * 'id ’ 

iS .- *J 



ro 

rs> 


TABLE II 

LIQUID STATE PHYSICAL PROPERTIES REFERENCE LIST (BOOKS) 


AUTHOR 

TITLE 

PROPERTIES 

MATERIALS 

McBride, Heimel, Ehlers, 
Gordon 

Thermodynamic Properties to 
6000® K for 210 Substances 
Involving the First 18 
Elements, NASA SP-3001 
(1963) 

heat capacity, enthalpy, 
entropy, free energy 

At ,AI 2 0 3 ,B,B 2 0 3 ,Be,Ba0 

li,lici,lif,li 2 o,lioh, 

Mg,MgF 2 ,Na,P,S, and SI 

Mantel 1 * Ed. 

Engineering Materials 

Hand*** ok 

(1938/ 

T , B.P., heat of fusion, 

ni 

density, vapor pressure, 
heat capacity, n, <e, o, 
surface tension 

Al ,Sb,BI ,Cd,Cs,Ga,Au, 
ln,Pb,LI ,Mg,Hg,K,Rb, 
Ag,Na,TI ,Sn,Zn 

Jackson, Ed. 

Liquid Metals Handbook 
(1935) 

isothermal compressibility, 
density, n, *, o, heat 
capacity, thermal dlffusl- 
vlty, Prandtl number 

Na 

BurdI 

SNAP Technology Handbook, 

NAA-SR-8617 

(1964) 

T . B.P. , density, n, «, o, 

ID 

heat capacity, surface 
tension, vapor pressure 

Hg,Na,NaK,K,Rb,LI 

Sharp 

Thermodynamic Functions 
for Water — — 

UCRL 7118 
(1963) 

P,V,T, entropy, enthalpy. 
Internal energy, Gibbs tree 
energy 

H 2 0 
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AUTHOR 

TITLE 

PROPERTIES 

MATERIALS 

Svehla 

Estimated Viscosities and 
Thermal Conductivities of 
Gases at High Temperatures, 
NASA T R-R-/32 
(1962) 

Lennarc -Jones constants, 
n , heat capacity 

o 

AI p AICI,AI0,AIS,BF 5 , 3 
B 0 ,H 2 0 ,Hg,H 9 Br 2 ,HgCl 2 , 3 
Hgl 2 ,L!,LIBr.LIO,U 2 , | 
LI^OjMg.MgCI ,MgF,Na, 




NeBr ,NaC 1 ,NeF ,Na 1 ,NeO , 
Na0H,0 2 ,P,Si,SICI,SIF, 




SIO,SIS,SI 2 ,X«,Zn 

Made lung 

Physics of lll-V Compounds 

atomic radii, T 

m 

B,C,N,AI,SI,P.Ga, 
Ge,As, 1 n, Sn, Sb,T 1 , Pb , 
Bl 

Glazov, Glagoleva, 
Chizhevskaya 

Liquid Semiconductors 
( 1969) 

electrical conduct Iv lty,o 

thermoelectric power, T*E.P. 
density, p 

kinematic viscosity, v 
magnetic susceptibility, x 
melting temperature, T M 
heat of fusion, Q. 
mobl 1 Ity, v 

SI ,Ge,Ta,AISb,Ga$b,lnSb, 
GaAs , 1 nAs ,ZnTe ,CdTe ,G« 2 T« ^ 

1 n 2 Tej,Mg 2 S 1 ,MQ 2 GdfMQ 2 Sn, 

MfljPb , GaTe , SnTe , PbTa , PbS* , 

Pb$,BI 2 $« 3 , 8 l 2 T« 3 # Sb 2 T * 5 

Hultgren, Orr, Anderson 
and Kelley 

Selected Values of Thermo- 
dynamic Properties of 
Metals and Al loys 
(1963) 

heat capacity, c 
enthalpy, Hy-H f p 
entropy, 

free energy, F^-H^ 

66 matats and samt -mats 1 s 
167 alloys 


atomic weight, A 
vapor pressure, *P 
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l\> 

fNJ 


TITLE 


PROPERTIES 


MATERIALS 


Gray, Ed. 


A IP Handbook 
' l 963) 


Isothermal compressibility, 

6 

surface tension, a 

thermal expansion coefficient, 

a 

thermal conductivity, « 


H^O , Hq , A I , Sb ,0 1 ,Cd ,Cu , Ge , 
Au,Pb,Hg,Pt,K,Ag,Na,S n ,Zn 
H 2 0,Hg,LI ,Na,Pb,Sb,Sn,Te,Zn 



Davydov 

German 1 urn 

heat capacity, C 



(1966) 

P 

melting temperature, entropy 

ro 



enthalpy: T m .S m .H m 

i 

00 



Lennard-Jones parameters, 
e4t> 

o 

molecular volume, 



1 


Beer, Ed. 


I 

! 


Liquid Metals - 
Chemistry and Physics 


thermal coefficient of AM metals and semt-metals, 

resistivity. Isothermal com- many alloys 
press! bl I Ity, thermal expan- 
sion coefficient, n, heat 
capacity, surface tension. 

Hall coefficient, resistivity, 
magnet I c suscept I b 1 1 1 ty 


r 4 
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2.2 (Continued) 


where L is a characteristic lengtn for the liquid system. To obtain 
values of n m in poise, the usual viscosity unit, the following con- 
versions are useo: 

1 »-* * &SHT - 


1 poise = 


dyne sec 

riP' 


The meaning of the magnetic viscosity is indicated by a dimensionless 
parameter, the Hartmann number, which is defined as 



where n is the molecular viscosity of the melt. Thus if h > 1, then 
the magnetic viscosity predominates, while for h < 1, magnetoviscosity 
is not as important. Figures 1, 2 and 3 are graphs of Hartmann number 
versus magnetic induction for various materials, values of L and, 
in the case of sodium, different values of temperature. These graphs 
show that, while liquid sodium is, as expected, most affected by a 
magnetic field, the III - V compounds I :Sb, InAs and GaAs do show a 
definite field-induced viscosity with h > 1 for fields > 40 gauss and 
L > 1 cm. A field of 40 gauss is quite modest. Thus magnetic vis- 
cosity dominates for most of the conditions shown in the figures. 

Convection Decay in a Magnetic Field 

The ilavier-Stokes equation describing bulk liquid motion in a magnetic 
field may be written 

p = Z + (J X g) (3) 

where ?. = -v [oii» + P - v-uj + n v 2 u contains all non-electromagnetic 
forces. From the expression for the current density, 3, in a moving 
medium and from appropriate vector identities, it can be shown that 
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FIGURE 1: HARTMANN NUMBER VS. MAGNETIC INDUCTION 

LIQUID SODIUM 
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2.2 (Continued) 

the component of magnetic force transverse to the magnetic field lines 
is 

(3 X B) t = - o B 2 u t . (4) 

If nr* considers only transverse components of equation (3) and defines 
a time constant, t, by 



equation (3) may be rewritten as 



(5) 


( 6 ) 


If it is assumed that the forces contained in change only slowly 
with time, (6) may be integrated to yield 


u t ■ li, * <c„ - i r ) e-V' 

t p t Opt 


(7) 


where c q is an integration constant equal to the free convection velocity 
at time = tg.the magnitude of which is on the order of 30 to 40 cm/sec. 
Figure 4 presents the qualitative decay curve for a situation corresponding 
to equation (7). Thus when a magnetic field is imposed on a melt in which 
convection currents are flowing, the flow velocity will be damped by the 
field according to equation (7), falling off after times much greater than 
r to a value 

“t {t y> ^ ih 

From equation (3), if there were no magnetic field. 



D256- 10024 


2.2 (Continued) 


*0 



p 


T 


t 



FIGURE 4: CONVECTION DECAY IN A MAGNETIC FIELD 


For liquids of density p = 6 ^p-, conductivity of o = 10 3 .Q * 1 cm* 1 and 
a magnetic induction cr 10 3 gauss, i is on the order of 6 seconds. 

Assuming that c Q ' 35 |^r, the velocity after three minutes will fall 
to about one cm/sec. 

Body Forces Due to Oscillating Magnetic Fields 

The qualitative aspects of oscillating magnetic fields have been dis- 
cussed in References 1 and 2. Since the body force in d ^r S on an element 
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2.2 ^Continued) 


of liquid is given by the cross product of the eddy current, 3g, 
and the magnetic induction, $, the magnitude of this "eddy current 
force" is 

% = |Ji l$l sin 0 (8) 

where e is the angle between the 3 and 5 vectors, f then assumes its 
maximum value when 9=5-, which is often the case experimentally, by 

6 (3) 

relating the eddy current to the field inducing it, one can show' ' that 

f p = jO | |B| = |WB = = ^2. (1+ X )H 2 (9) 

max 


where a is a characteristic length in the liquid, such as the radius 
of a cylindrical melt container. If one assumes that a is 5 cm, then 


2u 


= 0.05028 


dyne 
cm amp z 


Thus, the eddy current body force of equation (9) may be plotted as a 
function of field for different materials (different susceptibilities, 
x) as in Figure 5. In practice, since values of y for liquid metals 
and semiconductors are on tne order of 10 -6 cgs units, curves for the 
materials InSb, Ge, Na and Hg all fall on the same points to witnin 
0 . 001 %. 


To gain a more quantitative feeling for the magnitude of the eddy 
current body force, it is helpful to compare it to the magnetoviscous 
force arising from the magnetic viscosity discussed earlier and to 
the gravitational body force, pg. From equation (4), the magneto- 
viscous force may be written 

f m = o u Q 2 (1+x) 2 H 2 u (10) 
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2.2 (Continued) 

Thus by dividing equation (9) by (10) one obtain'. 



2 

o u 0 a(l+X/U 


( 11 ) 


where it must be remembered that fe max and are vectorial ly in 
different directions, and in fact are often orthogonal. Assuming 
again that a = 5cm, o = 10 3 iT^cm -1 and noting that 1+x : 1 and 
u 0 s 4n X 10" 9 n sec/cm, one may plot the ratio f 6max ^m f° r a ran 9 e 
of convection velocity values as in Figure 6. This graph shows that 
the predominance of the eddy current force increases for a diminishing 
or damped convective flow. Thus in using oscillating magnetic fields 
to control convection, or even to achieve stirring of the melt, one 
must be careful not to “overdamp" or set up eddy current-driven flows 
which are more vigorous Thar the gravity-driven convection at the 
beginning of the process. 


To determine just how the eddy current force compares to the force 
driving natural convection, equation (9) is divided by the gravi- 
tational body forv,e 


“max 



2tJ o (H *) 
a g p 


For a - 5cm, 1+ x = 1 , u Q * 4n X 10“ 2 and g = 980 , this 

equation becomes 


( 12 ) 



(5.13 X 10' 5 S-jr ) 
' amp^ crtr ' 


Hi 

p 


(13) 


which is plotted in Figure 7. It is obvious that f emiV dominates the 
liquid flow for field strengths above 500 |jjp- . Liquid sodium is the 
most responsive, of the materials considered, to the eddy current force, 
while the semiconductors InSb and Ge fall between sodium and mercury 
in their response to the field. Thus it is seen that oscillating magnetic 
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FIGURE 6: EDDY CURRENT/MAGNETOVISCOUS FORCE COMPARISON 
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2.2 (Continued) 

fields of even modest strengths (500 - 10,000 |®) produce forces 
which are quite powerful compared to other forces in melts. 

2.3 EXTERNAL FIELD EFFECTS DERIVED FROM THE FREE VOLUME MODEL 


The effects which magnetic and gravitational fields have on the 
temperature, volume, Gibbs free energy, diffusion coefficient and 
solidification rate of solidifying melts have been described qualita- 
tively in Reference 1. That analysis may be carried further by con- 
sidering the ratio of the solidification rate under changed field 
conditions, U', to the solidification rate under normal field conditions, 
U c - “Normal" field conditions implies zero magnetic field and 980 
gravity, while “changed" field conditions means either an applied 
magnetic field or microgravity conditions existing in an orbital 

(4) 

laboratory or other zero gravity simulator. Using Turnbull's expression' ' 
for the solidification rate (see equation 28 in Reference 1), and assuming 
that the fraction of lattice sites in a liquid-solid interface to which 
molecules can be attached remains constant, the field effect on 
solidification rate may be expressed as 



CT 

D 


l 

1 


giG'/kT* 

e AG/kT 


(14) 


The Free Volume Model expression for the diffusion coefficient, D, may be 
(51 

written yv* 

D = j 'i v* e” v^ (15) 

where 4 / is a parameter relating the distance traveled by a molecule 
between collisions to the average specific volume, v, of the liquid, 
and y is an overlap factor lying between j and 1 . v* is the critical 
value of free volume for the onset of diffusion, u is the gas kinetic 
velocity given by 
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and v,p is the free volume defined as 
v f = V - v 0 

where v Q is the molecular volume calculatec from the diameter and 
m is the molecular mass. 


( 17 ) 


Now if equations (16) and (17) are substituted into (15), one obtains 

j lll 


= $ T z e" v 


(1 - v^v )- 1 


(18) 


where 4> = v* is a constant for a given material and values of 

w* Jin ( 6 ) 

J -g- have been tabulated by Cohen and Turnbull' ' for some simple liquids. 


Since T and v are the only variables in equation (18), the temperature 
and specific volume under altered field conditions may be expressed 
as 

T' = T + aT and v' = v + av 


where aT and av are the temperature and specific volume changes produced 
by field changes. Writing the expression for D' in terms of T' and v' 
then dividing by equation (18) yieios the relation for the field effect 


on diffusion coefficient 

r- 


f(V,Av) 


(19) 


where 


f(v,Av) = (*£) (^) [(1 - 0)2 + M (1 . -° )r i 


Thus to determine the effects of magnetic fields or microgravity on 
U and D, the terms =!■ , — and the altered free energy, a6‘ , need 

V* I V 

to be calculated. 


Temperature, Volume and Free Energy Changes in Microgravity 



of isot 
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2.3 (Continued) 


change in the volume of a liquid due to a change in pressure, aP is 

= -BAP (20) 

where B is the isothermal compressibility of the liquid. If a container 
of liquid is envisioned as being moved from the Earth's surface into an 
orbiting laboratory where microgravity conditions prevail, thare will 
be a change in the hydrostatic pressure within the liquid of 

aP : - pZAg (21) 

where p is the liquid density, z is depth below the liquid free surface 
on Ea^'i and Ag is the change in gravity or acceleration field. Thus 
froir on (20) 

?- = BpZAg (22) 


From the definition of the thermal expansion coefficient, a, one can 

write 


aT av/v 
T ‘ «T 


(23) 


so the temperature change corresponding to the pressure-induced volume 
change may be calculated from (22) # 

(8 ) 

Since the Gibbs free energy of solidification is given' ' by 


AG = S fiT 
c 


.-where S c is the entropy of solidification and 
tooling (a negative number), externj 

aG' = aG + avaP 



im ' 20 ) 


a£ 
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with aP given by equation (21) for a microgravity case. W 1 M 1 these 
equations, the free energy term 

(1 - e aG */ kr ) / (1 - e iG/kT ) , 

the diffusion coefficient change D'/D, and the solidification rate 
change U^/U c may be calculated for any material if the characteristic 
parameters are available for the liquid state of the subject material. 

Table III lists values of the parameters required for calculations in 
both the microgravity case and in the magnetic field case for four 
representative materials. 

Temperature, Volume and Free Energy Changes in a Magnetic Field 

Equation (20) may be used to determine fractional volume changes in a 

(91 

liquid placed in a magnetic field if aP is now the magnetic pressure' ' 

aF> = " \ u o ^ H o ' xH2 ^ (27) 

where is the permeability of free space (.1257 dyne/amp 2 ), x is the 
magnetic susceptibility, H Q is the initial field strength before the 
liquid was placed in the field and H is the resulting steady-state of 
the field internal to the liquid. Thus 

r • \ 8»o [H| - x« 2 3 • (28) 

To estimate av/v without having to specify boundary conditions of a 
particular system, one may assume that for paramagnetic or diamagnetic 
liquids, H : H q . Since x is on the order of 10“ 7 for such liquids, 

^ j 8u 0 0 -x)H 2 = ^ ^‘ j 0 h2 • (29) 

An alternative method for calculating the volume change would be to use 
the theory of magnetostriction^^ . Unfortunately several parameters 
required in the equations of this theory are not available for liquids. 
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TABLE III PARAMETERS USED FOR CALCULATING FIELD EFFECTS* 




x s 

Q 

B 

a 

S c 

YV* 

V 

v o 

V 

T 

p * 

MATERIAL 

cm” ^ 

cm" 3 

erg/ 

cm 3 

2 

cm / 
dyne 

o “1 

K 

erg/ 

mole°K 


at?/ 

mole 

at?/ 

mole 

°K 

gm/ 

cm 3 

InSb 

-6 .816X1 0 -7 

- 7 . 612X1 0 - 7 

2.69X10 6 

3.15X10” 11 

.000341 

6.211X10 8 

0.1 

10.8 

18.2 

793 

6.43 

Reference 

13 

13 

11&14 

20 

20 

14 

E 

E 

39 


39 

Ge 

-4. 902X1 O -7 

-6.125X10' 7 

2.15X10 6 

1.19X10' 11 

.0001574 

2.443X10 8 

0.3 

6.8 

13.0 

1200 

5.57 

Reference 

13 

13 

11&14 

20 

20 

15 

E 

E 

39 


39 

Hg 

-2.466X10' 6 

-2.128X10” 6 

■ 1.52X10 5 

3.9X10" 12 

.00018 

9.8X10 7 

.0935 

8.5 

14.6 

224 

13.7 

Reference 

18 

18 

11&14 

16 

16 

14 

6 

E 

17 


17 

Na 

5. 859X10" 7 

4. 896X1 0" 7 

1.02X10 5 

1.916X10” 11 

.00028 

7.01X1 0 7 

.236 

17.3 

24.6 

361 

0.93 

Reference 

19 

18 

11&14 

17 

18 

14 

6 

E 

17 


17 


* H is assumed to be 10 5 oersteds (7.96X10 4 amp/cm) 

<5T is assumed to be -10°K (T * T m + 6T) 

m 

Z is assumed to be 1 cm 


E Indicates estimated values 
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In his paper describing magnetic field effects cn the dissolution and 
solidification rates of paramagnetic crystals in solution, Schieber^^ 
derived the following expression for the change in temperature of a 
solidifying system due to the application of a magnetic field, H 

r s (x * • x s ) W (30) 

where is the magnetic susceptibility per unit volume of the liquid, 
x $ is the susceptibility of the solid and Q is the latent heat of 
solidification. Thus the temperature at the interface of a solidifying 
material in a magnetic field can increase or decrease depending on 
whether the material is diamagnetic or paramagnetic and on the relative 
magnitudes of x £ and x s - 

The Gibbs free energy of a material in a magnetic field is given by 
Wood^^ as 

G' = G - vll * ft (31) 

where M is the magnetization, xH- Thus the change in 6 at the interface 
with a magnetic field applied will be 

aG' = aG - v(x£-x s )H 2 . (32) 

It is easy to show that if S C <5T > 4vQ, the free energy term 
(1 . e aG ' /kT ’ ) / (1 - e i6/kT ) 

is less than unity if x fc " x s is positive and becomes greater than one 
if x 4 - x s is negative. 

Calculation of Field Effects on Melts 

From equations (22), (23) ,(26) ,(29) ,(30) and (32) and the parameter values 
given in Table III, the various terms in equations (14) and (19) can be cal- 
culated and used to determine values for D'/D and U'/U for both the microgravity 

V \0 
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and magnetic field cases of altered external field conditions. Note 
that the references for the input data are indicated under the approp- 
riate values in Table III. Values of the variable parameters which 
were assumed for the examples presented here are given below the table, 

and estimated values of v^ and were determined from atomic and 

o v 

ionic radii shown on the Sargent-Welch Table of Periodic Properties of 
the Elements (1968). In such estimations a 50% composition for InSb 
is assumed. 

The values for U^./U c , D'/D and the free energy term which were computed 
from equations (14) and (19) are given in Table IV. Several comments should 
be made concerning the Free Volume Model and the results presented in 
Table IV. First, Turnbull's model is based on such simplifying assumptions 
that it could be expected to be most accurate for the simplest liquids. 

But when the model is compared with experimental data (see Table I, 

Reference 20), it is found to be more accurate for water and methanol 
than for such simple liquids as Argon ana Helium. In addition, Turnbull's 
purpose in developing the model was to describe transport and solidifi- 
cation in glass-forming liquids, but he also gets good agreement with 
experiment for organic liquids as well as liquid metals^. So there 
is little or no correspondence of the model to particular classes of 
liquids. On the other hand, because of the simplifying assumptions, 
the model is not overly accurate for any specific liquid^). But it 
can provide order-of-magnitude estimates for quantities of interest, 
and at least indicate the direction of change in these quantities 
under the influence of perturbing fields. This is the primary '-enefit 
of the model to Space Processing - that it does allow the determination 
of direct effects of magnetic and gravitational fields on diffusion 
coefficients and solidification rates. By "direct effect" is meant 
those effects derivable directly from the fields themselves, and not 
related to convective or magnetohydrodynamic effects. Such effects, 
shown in Table IV, can be expressed as percent changes if one notes 
that for numbers as close to unity as these, the percent change in 
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TABLE IV 

EXTERNAL FIELD EFFECTS ON REPRESENTATIVE MATERIALS 


MATERIAL 

MICROGRAVITY CASE 

MAGNETIC 

FIELD CASE 


FREE ENERGY TERM 

D'/D 

u'/u c 

FREE ENERGY TERM 

D'/D 

u '/u c 

InSb 

0.9999993 

1.0000005 

0.9999998 

0.99993 

1.00742 

1.007 

Ge 

0.99999966 

1.0000003 

0.99999996 

0,99987 

1,00624 

1.006 

Hg 

0.9999987 

1.0000007 

0.9999994 

1.00539 

0.99805 

1.003 

Na 

0.99999983 

1 .0000001 

0.99999993 

0.99771 

1,02139 

1,019 

1 
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a quantity X is approximately 100 (X'/X-l) where X' is the perturbed 
value. Expressed this way, the changes in both diffusion coefficients 
and solidification rate are in the ranges 10 -5 to 10* 6 ^ for the 
microgravity case and 0.3 to 1.9% for the magnetic field case. 

It must be noted that the changes in D and U c in the magnetic field 
case are strongly dependent on field strength. That is, for high 
enough fields, the value of D'/D for mercury could be greater than 
one, while for sufficiently low fields, the values of D'/D for InSb, 

Ge and Na could be less than one. but the data do show that gravity 
has a much smaller direct effect on D and U c than do magnetic fields. 
Thus this analysis supports the thesis that suppression of gravity- 
driven convection, which has an indirect effect on solidification 
parameters, is the primary benefit of materials processing in space. 

The data in Table IV also show that the free energy term is less 
sensitive to magnetic field changes than is the diffusion coefficient. 
This is due to the more direct dependence of the diffusion coefficient 
on field-induced volume change and to the fact that volume changes 
are much larger (by factors proportional to H 2 ) m magnetic fields 
than in microgravity where is on the order of 10 -8 . The increase 

in U upon the application of a magnetic field indicated in Table IV 

c (11) 
has been observed experimentally by Schieber' ' . 

2.4 THE RELATION OF EXTERNAL FIELDS TO SOLUTE DISTRIBUTION IN 
SOLIDS 

Previous sections of this report have dealt with calculations performed 
with theories developed during the period described by Reference 1. 

The present contract was primarily concerned with such calculations, 
but also provided for further theoretical development in the area of 
external field effects on solute distribution in solids. The question 
addressed by this portion of the work could be stated: "Is it possible 

to develop a theoretical model which will describe the effects an 
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external field has on the substructure of a material which solidifies 
in the field?" The best answer to this question at present is that 
a model consisting of four steps has been developed which formally 
relates the occurrence of microsegregation in a solid to the external 
field conditions which might obtain during solidification of a binary 
material. The steps are: 

1. Calculate the field distribution inside the melt due to the 
known external field as a function of susceptibility (magnetic 
field case) or density (gravitational field case). 

2. Find the force the internal field exerts on solute atoms. 

3. Determine the resulting solute atom velocity produced by the force. 

4. Use Sekerka's interface stability theory to relate the criteria 
for microsegregation to solute atom velocity. 

Although the model yields a formal relation between external fields 
and microsegregation, it is not really useful for performing practical 
calculations since each step contains operations which cannot be 
performed rigorously for most real solidification situations. However, 
it is instructive to consider the theoretical aspects of each of the 
above steps. 

Field Distributions in Melts 

The calculation of magnetic field distributions fi(r,t), in materials 

( 21 ) 

from a known field outside the material has long been understood' ' . 

basically, one solves Laplace's equation 

v 2 * = 0 (33) 

for the magnetic potential, $, where 

H = - v * (34) 


2-29 



0256-10024 


2.4 (Continued) 

in all regions of space subject to the boundary conditions of the problem. 
As in all field problems, the geometry of the situation determines the 
form of the boundary conditions and the solutions. Figure 8 shows 
the two geometries considered. These geometries were chosen because 
of their widespr^ao experimental usage. 



(A) 



(B) 


FIGURE 8: CRYSTAL-MELT GEOMETRIES CONSIDERED IN THE STUDY OF MAGNETIC 

FIELD DISTRIBUTIONS IN MELTS. (A) CORRESPONDS TO ZONE 
REFINING AND (B) TO CZOCHRALSKI CRYSTAL GROWTH. 
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The zone refining geometry can be considered an infinite cylinder, 
since the solutions of interest are those in a small region near the 
liquid-solid Interface. The boundary conditions in this simple case 
are 


(0,e) * finite 
ill, . 3 *o, 

W 1 ar 1 u 0 ar * 

(35a) 

(35b) 

a a 


♦l (a,e) = (a, e) 

(35c) 

♦q (r,e) — ► - H Q r cos e 

(35d) 


JP oo 


where $ 0 is the potential outside the cylindrical container, <j^ is the 

potential inside and is the initial field before the melt was intro 

0 (211 
duced into the field. The solutions obtained in the usual manner' ; 

are in cylindrical coordinates 


% - H q [cos e r-sin e e] + 

(^) fr H 0 [cos e r+sln e e] 

^1 = 2+7 H o ^ cos 0 r ‘ slR 6 

and in rectangular coordinates 

% = H 0 1 + (?+£) H o ?£ C( x2 -y 2 ) 1 + 2xy j], 
r 2 =x 2 +y 2 


a, * (£) « 0 i 


(36a) 

(37a) 

(36b) 

(37b) 


where i, j, r and 0 are unit vectors, with H Q = H Q i initially. 

',ie Czochralski geometry presents a much more difficult problem. First, 
because this case is essentially two cylinders of different radii 
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meeting end-to-end at the z=0 plane, there are more boundaries, and 
thus boundary conditions (see Appendix B) and solutions to be considered. 
Also, because the problem is not symmetrical on either side of the 
z=0 plane, the general solution to equation (33) i« a Fourier-Bessel series 

* * tVm< kr > * 6 km N m< k ' » X 

[C k e kz + D k e' kz ][E m cos me + F m sin me] (38) 

and this series unfortunately does not converge in the region near the 
z-G plane. Thus no rigorous analytical solution to this problem is 
possible, and one must choose either an approximate analytical solution 
or a numerical computer solution. An approach to obtaining an approximate 
solution is to convert the Czochralski geometry into the geometry of 
two infinite, concentric cylinders as shown in Figure 9. The solutions 
to the infinite concentric cylinder problem are easily obtained, then 
the solutions may be written above z=Q by setting = uq and below 
z=0 by setting ^ = Uj • 



* 2b * 


i 

X 2 


u l 

£ 



CRYSTAL 


y l 

1 

P 2 ! 

*1 

MELT | 

| 


— 2a *- 


FIGURE 9: CONCENTRIC CYLINDER APPROXIMATION TO 1ZCCHRALSKI GEOMETRY. 


2-32 



0256-10024 


2.4 (Continued) 

Thus, in cylindrical coordinates 

% * H 0 (cos e r-sin e e) + 

2 +^ 7? H 0 ( cos 0 6 0 ) (39a) 


whe^i x s X] and d = a, z < 0; x = x 2 and d = b, z > 0 


**1 = 2+x-j H 0 


(cos e r-sin e e) 


(40a) 


i? 


2 


l + xp 

H| 


0 


(cos e r-sin e e) 


In rectangular coordinates, the solutions are 

% = H 0 1 + 2+x H 0 F r + 2xyj] 


where x = and d = a for z < 0; x = x 2 and d = b for z > 0, r 2 = 
x 4 + y 2 and Aq indicates the new field outside the melt. 


fi l = 2+ Xl H 0 1 
l+x 2 

K = 2(^77=-) H Q i 


l ?+X, 


(41a) 


(39b) 


(40b) 

(41b) 


The trouble with these solutions is, naturally, that they do not 
satisfy the real boundary conditions (see Appendix B) in the liquid- 
solid interface region, which is the region of primary interest. An 
approach to obtaining j solution which does satisfy the boundary con- 
ditions for the field distribution in a Czochralski growth system 
utilizing computer capabilities is the finite element method. This 
method consists basically of four steps: 


1. break up the region of interest (interface region) into small elements 
as shown in Figure 10. 
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(A) 



FIGURE 10: FINITE ELE.SENT METHOD FOR SOLVING FIELD DISTRIBUTION 

PROBLEMS IN CZOCHRALSKI GEOMETRY 

(A) ELEMENT GRID IN INTERFACE REGION - SIDE VIEW 

(b) ELEMENT GRID IN INTERFACE REGION - VERTICAL VIEW 

(C) NODES OF ADJACENT ELEMENTS 


2. Solve Laplace's equation for <f> at each node for each element. 

3. Impose "compatibility" conditions on the nodes of element i so that 
solutions at each node of i satisfy the boundary conditions of each 
element adjoining i at that node, and therefore satisfy the boundary 
conditions of the problem. 

4. Program the computer to sum the solutions for each element (which 
satisfy boundary or compatibility conditions) over all elements to 
obtain the total solution of the problem. 

Techniques for actually carrying out the above steps are not yet per- 
fected since the application of finite element methods to field theory 
problems is relatively new. Thus the time and resources required to 
obtain a solution for the field distribution in a Czochralski crystal 
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growth system are far beyond the scope of the present study. However, 
it is possible to determine qualitatively the shape of the field in 
the interface region. Physically, it is known that magnetic induction 
(S) field lines are "pushed out" of materials possessing lower per- 
meabilities than their surroundings. It can also be seen from equations 
(40) and (41) that Bj > B£. Therefore the lines of B in the interface 
will be shaped qualitatively as in Figure 11 when xi > x 2 - Thus the 
qualitative field shape of Figure 11 coupled with the numerical answers 
from equations (39) - (41) constitute an approximate solution for the 
field distribution in the Czochralski growth case. 



FIGURE 11: B FIELD LINES IN THE INTERFACE REGION, CZOCHRALSKI 

GEOMETRY, X] > x 2 


Gravitational field distributions are even more difficult to calculate 
rigorously because they depend on the integrated matter distribution 
throughout the liquid-solia system as well as any external gravity 
field such as that of the Earth or the random acceleration field on a 
spacecraft. But oraer-of -magnitude estimates of gravitational forces 
on solute atoms may be found from F = mg^ where m is solute atom mass 
and g T is the total gravity field acting on it. On Earth g^ is the 
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sum of the Earth's mass field, its rotational acceleration field (these 
two sum to lg) and the gravity field of the materia, in which the atom 
is situated. 

For a melt in i cylindrical container, the gravity field due to the 

melt mass acting on a molecule near the surface of the melt may be 

( 22 ) 

estimated from Gauss's law v ' 

I / 9 c *ds = -4^ G / / / p dv (42) 

where G is the universal gravitation constant and p is the density of 
the melt. Assuming p and g c to be constant, equation (42) becomes 

V< 2r c + r c 0 ■ 
or 

9 C = -ttGp (r c r c + zK) (43) 

where r £ is the radius of the cylinder, z is its length and r c and K 
are the corresponding unit vectors. Now uGp = 1.26 X 10 -6 sec -2 , so 
for values of r £ and z on the order of 10 cm, the magnitude of g £ is 
roughly 10 -5 cm/sec 2 or 10 -8 g. Therefore, if r fi represents a unit 
vector lying along a radius vector of the earth, the total gravity field 
in a cylindrically shaped melt is 

g T = -g r e - 1 0~ 8 g (r £ + K) = -g r g . (44) 

In a space processing laboratory, the gravity level can be represented 
by 

g T = -g c + f (a) (45) 

where f(a) represents random spacecraft accelerations and vibrations 
which are on the order of 10" 3 g to 10 -6 g. 
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Resultant Internal Forces on Solute Atoms and Atom Velocity 

Once the magnetic field distribution in a melt is known the force on 
a given atom due to the field is^ 1 ^ 

? = ? (in! • B) 

mag i 


(46) 


where m! is the magnetic moment of the ith atom. Since in! is constant 
for a giver species of atom, equation (46) can be written 


f 


mag 


(m; 


7) 8 


or, in terms of scalar differences 

lb 


mag 


m 


lx 


(47) 


(48) 


The magnetic moment, m', may be calculated from the susceptibility of 
the material. Diamagnetic materials do not h&vp magnetic moments as 
such, but do exhibit behavior in a magnetic field which corresponds to 
an "effective moment" 

' ‘dia ■ T H 
0 


(49) 


where V is the volume per mole of diamagnetic liquid and N is Avogadro's 


number. For paramagnetic materials the moment is 

3 v Zr 0 - T 

"para = ^<r + T r ^ 

K 0 


(24) 


(50) 


where the second term on the right derives from the diamagnetic contri- 
bution to the paramagnetic susceptibility. If any solute atnm possesses 
a net electronic charge, it will also experience lorentz forces as it 
moves in the magnetic field of the melt. This situation is considered 

in Appendix C. 

Gravi tational forces, as have been indicated, are calculated from 


nig-, 


(51) 
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where m is the mass of a solute atom. Interatomic forces are foui.d 
from the negative gradient of the intermolecular potential, for 
the solute-solvent system in question. A value for the force may be 
calculated if one assumes that 2-body forces predominate and assumes 

(251 

a particular form for the potential. For the Lennard-Jones potential 

F mol = 24 r Kf) 12 - < 52 > 

where ( and c are the Lennard-Jones energy and distance parameters 
and r is the intermolecular separation. (The terms atom and molecule 
are used interchangeably since the simple models described here do not 
differentiate Detween molecules of different atomic structure.) Table V 
presents values of the forces calculated from equations 48 - 52 for 
several liquids. 

One must be careful, when comparing these forces, to remember that they 
are not static. That is, if a liquid is envisioned as a collection of 
molecules situated in cells formed by their nearest neighbors, each 
molecule is known to oscillate about the center of its cell with a 
restoring force such as equation (52). but since diffusion is known 
to occur, at some instant the force and cell structure must be altered 
enough to allow the molecule to escape its original cell and enter 
another. This alternation of the cell structure and restoring force, 
or cell potential, will occur randomly because, as must be recalled, 
each molecule making up the cell under discussion is also oscili ‘ing 
in its cell. Of course, even this picture is an oversimplifi cation 
because the description of liquid forces is one of the most difficult 
of the li-body problems. So although the data of Table V indicate that 
liquid intermolecular forces are mucn stronger than magnetic or 
gravitational forces, the statistical (and, in fact, unknown) nature 
of the true intermolecular force permits an external torce to act as a 
bias on the otherwise random motion of a solute molecule. 
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TABLE V + 

MICROSCOPIC FORCE COMPARISONS IN LIQUIDS* 


1 

Hg 

InSb 

- 

Na 

NH 3 

C°2 


- 

ergs 

1,0353 X in" 13 


1.898 X 10" 13 

7.707 X 10” 14 

2.694 X 10‘ 14 


a 

Ci.l 

2.827 X 10" 8 

-- 

3.567 X 10" 8 

2.9 X 10‘ 8 

3.941 X 10' 8 

CO 

r ~ 

cm 

3.59 X ID’ 8 

— _ 

4.27 X 10‘ 8 

4.3 X 10” 8 

5.02 X 10‘ 8 

cn 








UJ 
1 — 
LU 

m ~ 

gm 

3.33 X 10' 22 

1.97 X 10’ 22 

3.82 X 10’ 23 

2.82 X IQ' 23 

7.3 X IQ" 23 

§ 







et 

Xo ~ 

-3 

cm 

-2.466 X 10‘ 6 

-6.816 X 10‘ 7 

5.859 X ID' 7 


v mm 

CL. 

t. 







V— 

ID 

Cl 

h - 

w 

gm/mole 

200.6 

116.75 

23 

— 

-- 

21 

V 

cm^/mole 

14.6 

18,2 

2^ 6 

25 

40 


p - 

qm/cm 3 

13.7 

6,43 

0.93 

-- 

-- 


^mol 


8.6 X 10'° 

-- 

11.6 X 10" 6 

3.3 X 10' 6 

1.6 X 10‘ 6 

f 

t/0 oo 

F 


1.8 X 10" 21 

6.2 X 10' 22 

2.34 X 10“ 18 



LU Ll' 

mag 





O 21 
Q£ >- 
OQ 

F 


3.3 X 10' 19 

i io‘ 19 

3.74 X 10" 20 

2.76 X 10‘ 20 

7.15 X 10" 20 

Ll. 

9 






0 
ro 
c n 
o> 

1 

o 

o 

ro 

4 * 


* ff as*, .‘iied to be 600 ^ 

H assumed to be 60,000 oersteds 
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The next step in connecting external fields to solute distribution in 
the solid is to determine the solute molecule velocity, u, from the forces 
acting on the molecule. The practical method of relating velocity to 
force is through the mobility of the solute molecules 

u(r,t) = v ?(r,t) (53) 

where the term f includes all forces acting on the molecule except the 
random molecular diffusion forces, discussed above, which are included 
in the mobility, v. v depends on diffusion coefficient and temperature 
which are in turn related to the microscopic fluctuations, both random 
and oscillatory, of the molecules/^’^and therefore to the random inter- 
molecular forces. Thus, for the cases considered here, f would be a 
magnetic or gravitational force (as calculated from equations 46-51) 
while intermolecular forces, such as equation (52), would not appear 
explicitly in the expression for u. 

The Dependence of Microseqreqation on Solute Velocity 

Once the velocity of a solute molecule is known in the melt, how does 
this factor relate to solute distribution in the solid? The following 
theory, which discusses microsegregation in terms of the stability of 
the liquid-solid interface, provides a qualitative answer. 

Segregation may be divided into two parts: microsegregation and macro- 

segregation. Microsegregation includes short-range differences in 
concentration, such as those found between cells, dendrites and grains, 
while macrosegregation refers to long range variations in composition 
found along the length of an ingot and between the outside and the center 
of a casting. The difference in composition between the solid and 
liquid phases during solidification is responsible for the non-uniform 
distribution of solute in the final solid alloy. 
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The intercellular and interdendritic microsegregation, i.e., the 
segregation of solute to the cell boundary or inbetween dendritic 
branches, is the result of lateral diffusion of solute away from the 
tip of the growing projection. It follows that microsegregation occurs 
only when the planar solid-liquid interface becomes unstable. There- 
fore, to study the microsegregation and the substructure of the solute 
distribution, the morphological stability of the interface separating 
solid and liquid during solidification needs to be investigated. 

Interface instability in crystal growth was first discussed in connection 
with experiments reported by Rutter and Chalmers in 1953 on the 
unidirectional crystallization of dilute tin alloys in horizontal boats. 
They proposed that impurities rejected by the freezing solid can build 
up ahead of the advancing solid-liquid interface in such a manner that 
the equilibrium freezing temperature of the liquid adjacent to the 
solid-liquid interface is above its actual temperature. They postulated 
that if the "constitutional supercooling" exists, a protuberance on the 
interface would have a tendency to grow spontaneously, a smooth liquid- 
solid interface would be unstable and the observed growth forms would 
result. The conditions under which constitutional supercooling occurs 

e subsequently placed on a quantitative basis by Tiller, Rutter, 

( 29 ) 

Jacksc-i, and Chalmers' . The constitutional cooling principle deals 
w i l* r,c question of which state, solid or liquid, is thermodynamically 
wb:e in tne region of the liquid ahead of the interface and is not 
laV-o upon the dynamics of t.he whole system. 

A more elegant approach to the stability theory is to consider a small 
fluctuation on a planar interface and to determine under what conditions 
the fluctuation decays and under what conditions it grows. This is the 
per -.urbation theory of interface stability and has been developed by 
Wagner*''^, Mullins and Sekerka^ , Sekerka (33-36) an( j voronkev^^. 
The perturbation stability theory provides a description of the time 
evolution of a perturbed interface and the accompanying temperature 
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and concentration fields. This description is valid so long as the 
perturbations and their effects remain sufficiently small to be described 
by a linear theory. A time and distance scale of perturbed interface 
phenomena can therefore be ascertained^*^. 

There are two different but related calculational methods of studying 

the response of the interface to a perturbation. These are the time- 

independent and time-dependent stability theories. The time-independent 
(31-33) 

theory' ' involves the time implicitly but excludes time as an 

explicit variable in the description of the temperature and concentration 

fields. This "steady state” approximation gives valid results for 

sufficiently slowly growing perturbations. A more complex time- 

(34) 

dependent theory has been developed by Sekerka' ' employing time- 

dependent transport equations. The stability criterion of the time- 

dependent theory is very complicated and has not been reduced to a 

practical criterion. To an approximation, the necessary and sufficient 

condition for interface stability from the time-dependent theory has 

been found in agreement with that of the time-independent theory. 

Since the purpose of this study is to determine the effects of external 

forces on microsegregation of solute atoms during solidification, we 

. (32 33) 

have extended the time-independent theory of interface stability' ’ ' 

to include the movement of solute atoms resulting from the external 
forces which act on them. The stability criteria under different 
growth conditions are determined, and a method of estimating the size 
of the substructure is presented. 

Let us consider a unidirectional solidification system at constant 
solidification velocity, V as shown in Figure 12. The coordinate 
system attached to the planar interface moves with the constant 
velocity V with respect to the phases of either side. Suppose a 
sinusoidal ripple of infinitesimal amplitude, 6 , with a wavelength 
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FIGURE 12: INTERFACE GEOMETRY FOR SEKERKA'S THEORY 


x = 2 tt/ u satisfying the relation 

z = <fr(x,t) = 6 ( t ) sin w x 


(54) 


is introduced into the planar interface between liquid and solid. We 
wish to obtain an expression for 6 s d5/dt in order to see whether the 
ripple grows (6>0) or decays (i<0). This requires us to solve the 
following simultaneous equations, for the liquid. 


v 2 C + 1 (V-u) || = 0 (55) 

* * 4 s ■ ° < 56 > 


and for the solid 


7 2 T . , JL III = 0 

"ih 32 

subject to the boundary conditions 


T > = T i 


(57) 


(58) 


2-43 




D256- 10024 


2.4 (Continued) 


me + T m + T M rK 

<J> M M 

mC^ + 7 m - T M r6w 2 sin wx 


v(x,t) 


5 t K s <!P> - *1 <£> 1 

<p 

C (k-l) tD " S u -* 

4 > 9 


(59) 


(60) 


and to the conditions that for ( 6 ( much larger than any perturbation on 
the interface, T, T' and C shall take on the values of the solutions 
appropriate to the case of a planar interface. The notation is as 
follows: 


C 

0 

V 

u 

T 

T 1 


D 


th 


D ih 


k l ,k$ 

C L’ C s 


m 

L, 


Y 

Q 



= concentration of solute in the liquid 
= solute diffusivity of the liquid (diffusion coefficie.it) 

= average solidification velocity 
= solute velocity due to external force 
= temperature in the liquid 
= temperature in the solid 
K L 

= -=■ = thermal diffusivity of the liquid 
C L 

= K s = solid 
c s 

= thermal conductivity of the liquid and solid, respectively 

= specific heat per unit volume of the liquid and solid, respectively 

= slope of liquidus line on the phase diagram 
= absolute melting point of the pure solvent 

= y/Q,capillary constant 

= sol id- liquid surface free energy 

= latent heat of fusion of the solvent per unit volume 
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< = average curvature at a point of the interface 

v “ V + 3t = ’ nter ^ ace velocity 

k = distribution coefficient s ratio of the equilibrium concen- 

tration of solute on the solid side of the interface to that 
on the liquid side of the interface 

subscript $ denotes the quantity under consideration is measured at the 
interface. 


Let T and C be the temperature and concentration at the interface 

<t> $ 

respecti vely, we have 


T = T + a$ = T + a6 sin wx 
6 0 0 


C = C + bf = C + b6 sin wx 
6 0 0 


where T and C are the values for a flat interface and the second terms 
o o 

are the first-order corrections corresponding to the infinitesimal 

perturbation, a and b are coefficients to be determined. The solutions 

of equations (55), (56), and (57) satisfying the condition for large 

z and reducing to equations (61) on the interface z * $ are the following: 

G D . , i * 

C(x,z) - C Q = yT- [l-exp( : ^-z)]+ d ( b-G c ) sin wx e"“ (6i 


T(x,z) - T 0 = 


[l-exp(?p— z)] + <s(a-G) sin i»x e 


and on the solid side 


G 1 D 1 -r 1 z 

T'(x,z)-T 0 = y - - [1-exp (^ t—z)] + 6(a-G') sin w x e th 

th 


where G and G 1 are the thermal gradients at the unperturbed flat inter- 
face (a= 0) in the liquid and solid, respecti vely , G c is the concentration 
gradient in the liquid at the unperturbed flat interface, V';V-u is the 
solute velocity relative to that of the interface, and where 
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"* * ® * 1<®) 2 * « 2 i 7 , 

= + “ 2]7 

• - v 4 . r v 4 . 

“th ‘ ' 2 dT" [ 2DT" “ - 1 

th th 


i 

T 


(63a) 

(Ub) 

(63c) 


Substituting equations (62) into equations (59) and (60), we can 
determine a and b to the first order in 6 . From equations (62) and 
equation (60), after straightforward but tedious alget^a, we have the. 
expressio’ c or 6/6 which is sought 


6 

5 



mG c {.■>* - ^~)-[T^raj ? +^tG ' *- G)][^j* - jP(V - ^) ] 

2^ c ♦ | (V - jj)U* - § (V - jf)] 
k 


(64) 


as well as 

V = | (G'-G) * 1 (K $ G* - K l G) 



where 

- K s • K L 

K = _2__ — 1_ = average thermal conductivity of the system 
K L 

G = (— )G = aeneralized thermal gradient in the liquid 
K 


G* = (~)G‘ = generalized thermal gradient in the solid 
K 


(65) 

( 66 ) 


and p = 1 -k. 

Stability of the interface with respect to an undulation depends or. 
the sign of 6/6; a positive value of 6/6 for any u> means growth of 
some undulations and hence instability of the original flat interface 
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whereas a negative sign for all * means decay of all undulations and 
hence stability. In the p esent analysis, the conditions 0<k<l, m<0 and 
G‘ and 6>0 are assumed (because these hold true for the usual case). Also, 
only the condition of u>0, i.e., solute atom velocity in the same direction 
as the solidification velocity, is considered. There are two reasons why 
u<0 is not considered. Physically, u<0 means that the solute will pile 
up at the interface, producing a much more unstable ano thus more com- 
plicated situation than in the u>0 case. Mathematically, when u<0, the in- 
equality V-u > V- jj- (which will be utilized presently) does not hold. If 
this inequality is not valid, the denominator of equation (64) is not 
necessarily positive, and the stability analysis becomes too complicated 
to be accomplished within the time frame of the present study. But for 
the u>0 condition, three different cases are analyzed below. 


(a) 


V - - = 0 

P u DG 

From equation (66), we have V * — = - which weans G c -*0 and 

* ' n 



_D_ 

PC„ 


as V - - 
P 


0. Under this condition, equation (64) 


becomes 



mpC 


M“ 


{ (G'+G)]** 


. i 

K m * )C 0 


Everything on the right side of the above equation is positive except 
m which is less than zero. Therefore, we have found, when V - ~ = 0, 
that 
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which means that the flat interface is always stable under this 
condition. The physical meaning of V - jj- * 0 can be explicitly seen 
from the following relations. As V - u/p*0, or G c -*0, from equation (62) 
we have C(x,®) = C Q . This means that the concentration of solute at 
large distances from the interface is equal to the concentration at 
the flat interface, or, there is no pile-up of solute at the interface 
This condition can be realized by applying a force (such as a magnetic, 
acceleration or other external force as in equation 53) on the solute 
atoms in the direction of solidification with a magnitude of 

F * ar ey 5 ' - V 51 • (67) 


Hence, an external force acting on solute atoms can be expressed in terms 
of measurable material properties. 

(b) V - > 0 

Under this condition, the following relations may easily be found. 


and 


V‘ = V- u>V-^->0 since p s 1-k < 1; 

DG 

V - - = - > 0 hence G„ < 0 and mG^ > 0; 

P pCi c c 

o 

* .. )LL = V.~M > 2. ( '_ u \ > £ (u «. H.) 

“* D D D ' ; 0 ' p ; 


With these relations, all terms in the denominator of equation (64) 
are positive. The sign of 6/6, or the stability criterion, is solely 
determined by the sign of the numerator of equation (64) or the sign 
Of N(w) which is, from equation (64), 


kV 

NU) = " G c 0 - . -T7r~?r ] - T M r “ 2 - ? 

“o (v • p) 


( 68 ) 


= f( w ) - ^ (G'+G) + mG c 
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2.4 (Continued) 


where 


f ( uj ) r - T u ra) 2 - 2kmG 

n C 


v 

)T 


[1 


+ /2D 

+ (yr uj 


) 2 ] 


i/2 


- l+2k 


V 

\T 


(69) 


The function f(u) has the following properties: 

(i) f(w) < 0. 

( i i ) f(0) = -mG c and f(°°l = -<*> 

(iii) f(«) either decreases monotonically as u increases from zero 
or it has only one maximum for some positive value of uj. 


Properties (i) and (ii) follow immediately from the expression for f(w), 
equation (69). Property (iii) has been proved in a similar way as in 
Reference (32). The negativity of f(m) means that it always favors 
stability; evidently it poses the least barrier to instability at the 
frequency for which it attains its largest value (value of smallest 
magnitude). 


For the case in which f(<o) decreases monatonically as u increases from 
zero so that jmax f(w)j = | f (0) | = mG c> the sign of N(u>) is always 
negative. Therefore, the^e can be no instability under this condition. 
This absolute stability situation occurs when df/du <_ 0 for w -*■ 0. 

For small u, equation (69) may be expanded in a Taylor's series 

mG D 2 

fU) = -mG c + (^yry- - T^r )u 2 + .... 
then the absolute stability (df/du <_ 0 as u ■* 0) occurs when 
mG D 2 

WT * V - 0 


or 

(kyr) T f .:V 2 
mG c i) 7 


(70) 
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2.4 (Continued) 

Since, from equation (66) 

0G C • -PC 0 (V - £> 
and from equation (62) 

C(x,») = = C 0 k Jr 

the condition of absolute stability (equation 70) in terms of measurable 
quantities is 



which means that the capillary effect dominates the solute effect. 

(291 

The old constitutional supercooling criterion' ' is 

-G+mG c < 0 stable (72a) 

-G+mG c > 0 unstable (72b) 

Comparing 72b with the criterion of absolute stability (equation 70), 
ve may have a condition that 
kT M rV'V 

Z<^ c < -§2 (73) 


That is, the constitutional supercooling criterion for instability and 
the condition for absolute stability could be satisfied simultaneously. 
Hence, according to this theory, the constitutional supercooling is 
not a sufficient condition to determine the interface stability. 
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2.4 (Continued) 


For the condition that 

u nyv» 

p T^TopT 

or 

a (k yr) T M rV ' 2 k2 (\p") T n r v' , 

^ ra"; ■ 


(74) 


the stability criterion is very complicated. The stability conditions 
will be discussed below using an analytical method similar to that 
described by Sekerka ^ 33 ^ From equation (68), the stability criterion is 


N(u>) < 0 


or 


G' + G 


> 


1 + 


thl 

mG c 


From equation (69), we have 

,r i*) 2 


. = ii 


2 (kjr) 


(Sb 

C 


mG 


n + c 


2D 


TJ2 \ 

)] ' - 1 + 2 H 


(75) 


(76) 


A new function Y(y) is now defined by 
Y(y) = - lisil . 


mG. 


(1 + ey) 


2k’ 

TJT 


1 + 2k 1 


where 


and 


V 2 
> = ar“ 

k' = k L 

K -Kyi 


4k' 4mG c° 2 

a " yipT- 


(77) 
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2.4 (Continued) 

The criterion for a stable interface, from equation (75), becomes 


G’ + G 
2mG~ 


> 1 - Y(y) 


( 78 ) 


From property (iii) of f(w), the Y(y) may have a minimum at a certain 
value of y. Therefore, the system which is stable for all frequencies 
has to satisfy the condition 

G* + G 


O P ^ 

zmb c 

where S i l-[min. Y(y)] is defined as the stability function. 
To find the value of y at which Y(y) has a minimum, we set 


dY(y) , 
dy y = y 


= 0 


m 


From equation (77), we have 

. l/2 


0 ♦ ey m ) 


ek' 


[(1 + ay m ) 7 ' - 1+ 2k’] 2 


Since u must be real, we have 

(1 + ey ) J? - 1 + 2k' > 0 


(79) 


(80) 


and 

i/4 

(1 + %) > 0 

By taking the square root of both sides of equation (80), there follows 

r- + (2k' - l)r - (2k'/A 1/2 ) = 0 (81) 

where r (1 + ey W 4 . Since 2k'/A l/2 > 0 and the quadratic term is 
m 

missing, equation (81) has one and only one positive and real root of r. 

It is this positive root which corresponds to positive u. Since 
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2.4 (Continued) 

y m = ( r4_1 ) Ay4k ’ 

it follows from the definition of S that 
S(A.k') = l-[min.Y(y)] = 1 


(82) 


(83) 


According to equation (79), stability obtains when 


G' + G 

2mG > s ( A ’ k ') 
c 

Using equations (65) and (66), the stability condition, equation (84), 
can be written in terms of experimentally measurable quantities as 


2K, 


V k l 



D 



S(A,k* ) 


(84) 


(85) 


If all the material constants, K^, K $ , Q, h, y, D and m are specified 
by a choice of base material and solute, the lef f hand side of equation 
(85), defined as a test function, is a function of four variables of 
operation 


J 



( 86 ) 


Suppressing the explicit appearance of the material constants in the 
stability function, we call 

s ( F’ r ] = s ( A > k ') (87) 

oo 

Then for stability, we need the test function to be greater than the 
stalility function, i.e., 

J (£, }r, CJ > S (£r, ^ ) . (88) 

au 

By comparing the stability function, S(A,k') (equation 83), and equation 
(81) (which we have calculated including the solute atoms moving with 
velocity u) with Sekerka's stability function, (equations 11 and 14 in 
Reference 33), the functional form of both functions is seen to be 
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2.4 (Continued) 

identical. Sekerka^ 33 ^ has calculated his S(A,k) as a function of A 
in the range 0 jA< 1 for 16 values of k. Hi. results are plotted both in 
linear and semi logarithm scales as shown in Figure 2 and Figure 3 of 
Reference 33. His results can be directly applied to our results by 
changing his k and V to k'(;ky-r) and V ' (= V-u ) . 

(c) V - jj- < 0 

This condition implies, from equation (66), that the solute in the 
liquid is depleted at the interface (G c > 0) even for k<l. It can 
be easily seen that, from equation (64), 6/6 is always less than zero 
if G' + G is positive. Therefore, the planar interface is a stab’e 
condition during the solidification, and rnicrosegregation cannot occur. 

It is possible to calculate and plot 5/6 as a function of frequency 

(w=2tt/\) for any particular experimental case from equation (64). A 

schematic graph is shown in Figure 13 for a cell growth experiment. 

It can be seen that within the frequency range, 6/6 is 

positive and has a maximum value at w . That :< =»t oj=uj , the 

mm 

perturbation 6 has the fastest growth and the * cell size should 

be dominated by this perturbation. Therefore, >;e can estimate the 

order of magnitude of a cell size as .• =2-/j . 

3 m m 



FIGURE 13: 6/6 VS. FOR UNSTABLE INTERFACE CONDITIONS 
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2.4 (Continued) 

(32) 

As mentioned in Mullins and Sekerka's paper' ’ , there is no detailed agree- 
ment between the theoretically predicted relationship and the experimentally 
observed relationship of the cell size to the growth parameters. They 
gave two reasons for this failure: (1) time- independent or steady-state 

diffusion and temperature fields have been assumed in the calculation, 
and (2) the cellular structure which develops in actual experiments 

does not have a small amplitude. The first deficiency of the theory 

(34) 

was removed by Sekerka in 1967' ' . He formulated a time-dependent 

stability theory by using the same small perturbation technique. The 

results of that theory are very complicated. However, the cell size 

(35) 

predicted in a simplified version of the time-dependent theory 
does still not agree with the observed structure. The main problem is 
the second failure (above). To remove this obstacle, a finite amplitude 
theory rather than the small perturbation technique is required. Then 
the field equations and the boundary conditions all become nonlinear. 

To date, to our best knowledge, no solution to the nonlinear problem 
has been achieved. 

Thus our theory, which describes an external force field acting on 

solute atoms, follows a time-independent analysis similar to that of 
(32' 

Mullins and Sekerka' ( We can use the results to estimate the final 
cell size only to an order of magnitude. We could extend our present 
theory to a time-dependent theory, but even though such a theory is 
very complicated, it does not improve the prediction of the final solid 
substructures. 
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3. "Qualitative Effects of Osci ~ic Fields on 

Crystal Melts", Journal of Crys.-'. u. wt '0, 310, 
November 1973. 

4. "Thermodynamic Properties Derived From the Free Volume 
Model of Lioi'ids", Metallurgical Transactions, 5_, 643, 
March 1974. 

5. Further Analysis of Field Effects on Liquids and 
Solidification , Boeing Document D256-10024, July 1974. 

6. "External Field Effects on Solidification: Ma' oscopic 

and Microscopic Models", Proceedings of the 1974 AIAA/ASME 
Thermophysic; and Heat Transfer Conference, July 1974. 

7. “External Field Effects on Diffusion and solidification 
Derived from the F.ee Volume Model", Submitted to Journal 
of Applied Physics. 
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3. "Some Aspects of Present and Future Research in Space" 
(Invited Papei,, Annual Meeting of SESAPS*, November 
1973. 

4. “External Field Effects on Solidification: Macroscopic 

and Microscopic Models" AIAA/ASME Thermophysics and Heat 
Transfer Conference, July 1974. 
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APPENDIX B 

BOUNDARY CONDITIONS ON MAGNETIC FIELDS 
FOR THE CASE OF CZOCHRALSKI CRYSTAL GROWTH GEOMETRY 

c rom the Maxwell equation for the divergence of the magnetic induction 
field, v-B = 0, it can be shown^^ that the normal component of B must 
be continuous across each boundary, or 

(B 2 - B^ • n = 0 (B-l) 

where n is a unit vector normal to the surface. Since 

B = -J = v ; (B-2) 

this condition implies that 

-1 • =1 ■ '2 f '2 < B ' 3 > 

For this condition to hold, it is obvious that the potential, c, must 
itself be continuous across a boundary 

= : 2 (B-4) 

and to insure physically reasonable results, we must impose the con- 
ditions that : be finite everywhere and that at distances far from the 
boundaries of interest, the field approach the value the field would 
have if there were no material boundaries present. 

The folloving boundary conditions are based on the above principles, 
and the assumption that the cylindrical geometry is uivided into 
four regions as shown in Figure B-l , even though regions 3 and 4 are 
physically identical. 

(0,e,-«), : 2 (0,n,+~) = finite 
:{r,-;,z) ► - H Q r cos r 

f* -► Jt 


(«-5) 

(B-6) 


B-l 
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^ 3 (a,e,z) 

= $-j (a, 9 »z) 

(B-7) 

tii i 

O 3r ' 

a 

= . !!ii 
**1 > 

a 

(B-8) 

♦ 4 (b,0,z) 

~ i 2 (b,e,z) 

(B-9) 

u O 3r 

3^2 

y 2 *ri b 

(B-10) 

r ,s ,0) 

= $-|(r,0,O), b<r<c 

(B-l 1 ) 

U ° VT'o 

3$ 1 

* "1 arl Q - b<r<a 

(B-l 2) 



FIGURE B-l: CZOCHRALSKI CEOMETRY FOR MAGNETIC POTENTIALS 
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APPENDIX C 

MAGNETIC FIELD FORCES ON MOVING 
CHARGED MOLECULES IN A MELT 


A molecule which possesses a net electric charge, q, and which is in 

motion relative to magnetic field lines in a melt will experience 

Lorentz forces just as will any charged particle moving in a magnetic 

field. The evidence that particles in melts of semiconductor materials 

(39) 

are charged is stated by Glazov, et.al. , who suggest that in many 
such melts, q « 4e where e is the electron charge. Thus the Lorentz 
force on the molecule is 


= q u X 6 


(C-l) 


where 


u 




(C-2) 


with Uq being an scillatory velocity, r^, about some center of vibration, 
r, which moves through the melt with the translational velocity, u t =r t - 
The situation is represented schematically in Figure C-l. 


Center of 
Vibration 



Molecule 

Position 


3 


FIGURE C-l : ITINERANT UiCILLATOR GEOMETRY 


The velocities u Q and u^ may be obtained formally from the Itinerant 
Oscillator model ^). This model is based on the solution of the two 
coupled acceleration equations 


+ n o f 0 + a 6 (r o 


? t ) = A 


(C-3) 
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\ + t K = 6 (C ‘ 4) 

where ^ and ^ are friction constants, ^ is the oscillation frequency 
which is related to the intermolecular forces acting on the molecule, 

A is a randomly fluctuating driving acceleration also related to inter- 
molecular forces (and thus distances, a.) and the term 3 contains the 

J 

random Brownian motion-type acceleration and, after the magnetic field 
has been applied for some time greater than the initial transient 
(or build-up) period, also contains accelerations arising from the 
Lorentz force, equation C-l. Accelerations derived from the oscillatory 
motions should have a net time averaged value of zero. B must also 
include accelerations due to electrostatic forces between charged 
molecules. 

A more detailed analysis of magnetic forces on moving charged molecules 
in a melt is beyond the scope of this study, since such an analysis 
simultaneously involves two extremely difficult problems. Although 
the Lorentz force itself is well understood, actually calculating 
trajectories and velocities for particles which the force affects and 
which are moving through a medium other than vacuum is not so simple. 

The description of molecular motion (diffusion), is reasonably well 
understood on a statistical or macroscopic basis (see sections 2.3 
and 2.4), but on a microscopic basis, the motion of an individual 
molecule in the absence of external forces depends on random processes 
which are difficult to describe mathematically. 

Thus even though the Itinerant Oscillator model may provide a formal 
relation between Lorentz forces and charged particle motion in a melt, 
much more theoretical effort is required in order to determine the 
real nature of in melts of semiconductors. 
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